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Abstract. Recent discovery of Dsj particles, which are considered to be a great trouble by exper- 
imentalists as well as theorists, had already been resolved by our potential model proposed some 
time ago by two of us (T.M. and T.M.), in which the Hamiltonian and wave functions are expanded 
in l/(heavy quark mass) respecting heavy quark symmetry. 

Using our model, we explain how narrow states like D^j can be realized, predict properties of 
0+ and 1+ states of B and heavy mesons, and interpret how global SU{3) symmetry seems to be 
recovered for these 0+ and 1+ heavy mesons. 



INTRODUCTION 

Recent discovery of narrow meson states D5y(2317) by BaBar Jll], Dsj{2460) by CLEO 
i] and the following confirmation of both states by Belle [3^ has driven many theorists 
to explain these states since the former study of these states using a semi-relativistic 
potential model 0, ^ seems to fail to reproduce these mass values. This triggered a 
series of study on spectroscopy of heavy mesons again. To understand these states, an 
interesting explanation is proposed by Bardeen, Eichten, Hill and others [6, 7] who used 
an effective Lagrangian with heavy quark symmetry combined with chiral symmetry of 
light quarks but they cannot reproduce absolute mass values of these particles. 

This problem of Dsj has been considered by many experimentalists as well as theorists 
to be a great trouble until recently. It turns out, however, that actually we have already 
predicted these masses within one percent accuracy at the first order of purturbation in 
1/mg with mg heavy quark mass some time agoL8] by using the same potential model 
as in 10, m. Since the data adopted at the time of publicationJi] is obsolete, we are 
now refining calculations up to the second order. [9] The main difference between our 
treatment of the potential model and others is that we have taken into account negative 
energy states of a heavy quark in a bound state while others do not. Namely mass is 
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FIGURE 1. Procedure how the degeneracy is resolved in our model. 



expressed as an eigenvalue of a Hamiltonian of the heavy-light system, in which a bound 
state consists of a heavy quark and a light anti-quark and negative energy states appear 
in the intermediate states of a heavy quark in calculating an energy eigenvalue, that 
is missed in the former study. These negative-energy state contribution is not so small 
which is the main reason why people of [4, 5] could not get a right answer. 

Figure 1 . shows how our model generates masses step by step. At first chiral symmetry 
is intact with a limit of light-quark mass and scalar potential going to zero. At this stage 
no heavy quark mass correction is taken into account except for the lowest contribution, 
niQ, to a heavy meson mass, i.e., heavy quark mass is finite even at this moment. 
Next chiral symmetry is broken by introducing light-quark mass and scalar potential. 
There is still heavy quark symmetry remained so that some states with partial angular 
momentum, L = 0, 1, between light quark and heavy quark are still degenerate. Finally 
all degeneracies are resolved by including 1/mg corrections. The similar scenario is 
adopted in |0] except for the first stage. 



OUR POTENTIAL MODEL 
Dsj Particles 

Our formulation using the Cornell potential is to expand Hamiltonian, energy, and 
wave function in terms of 1 / mg and sets coupled euations order by order. The non- 
trivial differential equation is obtained in the zeroth order, which gives orthogonal set 
of eigenfunctions, and quantum mechanical perturbative corrections to energy and wave 
functions in higher orders are formulated. Consistency equation is given by 

4- //o + //i + ---)(^^o + ^n + ■■■)= (£o + ^[ + ---)(^« + ^^i + ---)(l) 

where integers of subscripts and superscripts denote order in l/mg, H = //fwt — f^Q, 
and the Foldy-Wouthuysen-Tani transformation is operated on a heavy quark and the 
Hamiltonian. [8J 

In Tables 1 and 2, stands for total spin and parity, Mq lowest degenerate mass, 
c\/Mq the first order correction, Mcaic calculated value of mass, and Mobs observed mass. 



TABLE 1. D meson mass spectra (units are in MeV). 



State f^'+'L?) 


k 


jP 




ri / Mr\ 






% 


-1 


0- 


1782 


0.477 xlO-i 


1867 


1867 




-1 


l- 




1.272 xlO"' 


2009 


2008 


'Po 


1 


0+ 


2073 


1.024 X 10-' 


2285 


2308 


"M" 


1 


1+ 




1.687 xlO-' 


2423 


2427 


"iPi" 


-2 


1+ 


2123 


1.412 xlO-' 


2423 


2420 


'P2 


-2 


2+ 




1.616 xlO-' 


2466 


2460 




2 


1- 


2329 


1.863 xlO-i 


2763 




"3£)2" 


2 


2- 




2.021 xlO-' 


2800 





TABLE 2. Dj meson mass spectra (units are in MeV). 



State C'+^Lj) 


k 


7^ 


Mq 


ci/Mo 




Mobs 


% 


-1 


0- 


1902 


0.351 xlO-' 


1969 


1969 




-1 


1- 




1.101 xlO-' 


2112 


2112 


'Po 


1 


0+ 


2102 


1.088 xlO-i 


2330 


2317 




1 


1 + 




1.762 xlO-' 


2472 


2460 




-2 


1 + 


2242 


1.269 xlO-' 


2527 


2535 




-2 


2+ 




1.462 xlO-' 


2570 


2572 




2 


1- 


2349 


2.006 xlO-' 


2821 




"3£)2" 


2 


2- 




2.170 xlO-' 


2859 





In Tables 1 and 2, we give our refined first order calculations both of D and Ds heavy 
mesons with various quantum numbers. Especially 0+ and 1+ states should be compared 
with experiments, i.e., D*0(2308), d'0(2427),|110|] A.y(2317), and Ay(2460).|Iilll The 
calculated masses, Mcaic, are within one percent of accuracy compared with the observed 
masses, M^/^. In the calculations, we have used values of parameters listed in Table 3. 



0+ and 1+ States of and B, 

In Tables 4 and 5, we give our refined calculations both of B and heavy mesons with 
various quantum numbers. For 5 and B^ mesons, unfortunately there are only a few data. 
We predict the mass of several excited states which have not yet been observed. Among 
them, calculated masses of 0+ and 1+ states for Bs mesons are below 5^/5*^ thresholds 
the same as the case of Ds mesons. Therefore, their decay modes to 5(0^/1^) +K/K* 
are kinematically forbidden, and the dominant modes are the pionic decay: 

Bs{0+)^Bs{0-) + n, 5,(l+)^5,(r) + ;r. (2) 

The decay widths of these states are expected to be narrow as for those of mesons, 
since these decay modes violate the isospin invariance. 



TABLE 3. Most optimal values of parameters. 



Parameters ttj 


a (GeV- 




/? (GeV) m„,,/ (GeV) m, (GeV) 


Mc (GeV) 


nih (GeV) 


0.2620 


1.937 




0.07091 


0.00803 0.09270 


1.040 


4.513 




TABLE 4. B meson mass spectra (units are in MeV). 




State C'+^Lj) 


k 




Mo 


ci/Mo 




Mobs 


% 


-1 


0- 


5255 


0.373 xlO-2 


5275 


5279 




-1 


r 




0.994 X 10-2 


5307 


5325 




1 


0^ 


5546 


0.882 X 10-2 


5595 






1 


l^ 




1.453 xlO-2 


5627 






-2 


l^ 


5596 


1.235 xlO-2 


5665 




^2 


-2 


2^ 




1.413 xlO-2 


5675 






2 


r 


5802 


1.723 xlO-2 


5902 




"3£)2" 


2 


2- 




1.870 X 10-2 


5911 





The similar decay modes are predicted for 5(0+) and 5(1+) mesons. 

B{0+)^B{0-) + 7t, B{l+)-^B{\-) + 7t, (3) 

but with broad decay widths, the same as ^(0+) and -0(1+) since isospin invariance is 
not broken and there is no threshold. 

These higher states of B and B^ mesons might be observed in Tevatron/LHC experi- 
ments in near future by analyzing above decay modes. We hope that our framework is 
confirmed in the forthcoming experiments. 

Recovery of Global SU{3) for 0+ and 1+ Heavy Mesons 

When one carefully looks at Tables 1 and 2, one notices that it looks that the global 
SU (3) flavor symmetry is recovered when = 0+ and 1+ for D(^s) mesons compared 
with 0^ and l^. The magnitude of SU{3) symmetry breakdown for D(^,-j(0^) and 
Z)(j-j(l-) mesons is given by (units in MeV), 

M{D{cu)) ^ M{D{cd)) ^ M{D{cs)) - 100, (4) 

while the symmetry seems to be recovered for £>(0+) and £>(!+) mesons, 

M(D{cu)) ^ M{D{cd)) ^ M{D{cs)) - (9 ~ 33). (5) 

For instance, mass difference among members of a multiplet (D(0^),Ds(0^)) seems 
to be larger compared with that among those of (D(0+), A(0+)). These equations are 
both for observed and calculated masses.ljsl ||9|] What causes this recovery of SU{3) 
symmentry? Is this just an accidental or is there a deep meaning for this? We study this 



TABLE 5. Bs meson mass spectra (units are in MeV). 



State C'+^Lj) 


k 




Mq 


ci/Mo 




Mobs 




-1 


0- 


5375 


0.287 


5391 


5369 




-1 


l- 




0.898 xlO-2 


5424 






1 


0+ 


5575 


0.945 xlO-2 


5627 




"M" 


1 


1+ 




1.531 xlO-2 


5660 




"iPi" 


-2 


1+ 


5715 


1.147 xlO-2 


5781 




'P2 


-2 


2+ 




1.322 X 10-2 


5791 






2 


1- 


5822 


1.866 xlO-2 


5931 




"3£)2" 


2 


2- 




2.018 xlO-2 


5940 





phenomenon in details lEl] and have found that this is a general phenomenon occuring 
for heavy mesons, (Qq), as well as heavy baryons like (QQq). 

This is explained by drawing a figure of mass gap between degenerate states, 
Mo{0-) = Mo{l-) (L = 0) and Mo(0+) = Mo(l+) (L = 1). These degenerate masses 
are written as Mq = mQ + Co{mq), i.e., a function of light quark mass niq. Hence defining 

AMK) =Mdo{0+) -Mz)o(O-), AM(m,) =Mz),o(0+) -Mo,o(0-), 

these mass gaps have the same values for D and B mesons. This mass gap function has 
monotonous decreasing tendency so that AM{mu) > AM{ms) for m„ < m^. 

AM(m,) - AM(m„) = -91 MeV, (6) 

which almost cancels mass difference between D{cu) {D{cd)) and D{cs) given by 
Eq.(4). Equation (6) is the origin of recovery of SU (3) invariance, which does not depend 
on whether the light-heavy system is Qq or QQq when q = u,d^ s. That is, this is a general 
phenomenon not peculiar to £>(0^) and D(l^). Relation between this phenomenon and 
constituent quark mass will be discussed elsewhere. Hllh 
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